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I.  INTRODUCTION 

Consider  a  perfectly  conducting  surface  S  placed  in  free  space  and 
illuminated  by  an  incident  electromagnetic  field.  The  problem  is  to  find 
the  surface  density  *J  of  electric  current  induced  on  S.  A  solution  for  £ 
can  be  obtained  by  writing  the  electric  field  integral  equation  on  S  and 
then  numerically  solving  this  equation  by  means  of  the  method  of  moments. 
This  solution  is  called  an  E-field  solution.  E-field  solutions  have  been 
applied  to  a  rectangular  bent  plate  [1],  a  surface  of  revolution  [ 1 ] — [ 5] 
and  a  surface  of  arbitrary  shape  [6]. 

These  E-field  solutions  are  reasonably  accurate  in  the  resonance  region. 
The  resonance  region  is  the  range  of  frequencies  for  which  the  maximum  di¬ 
mension  of  the  surface  S  is  of  the  order  of  a  wavelength.  Unfortunately, 
an  erroneous  change  in  the  slope  of  the  radially  directed  electric  current 
density  near  the  center  of  a  conducting  circular  disk  of  radius  0.02A  was 
reported  in  [5].  Here,  X  is  the  wavelength.  For  the  dual  problem  of  the 
circular  aperture  of  radius  0.02A  in  an  infinite  conducting  plane,  an  er¬ 
roneous  change  in  the  slope  of  the  radially  directed  magnetic  current 
density  appears  near  the  center  of  the  aperture  in  [7,  Fig.  6d] .  In  the 
E-field  solution  of  [5],  the  erroneous  change  in  the  slope  of  the  radially 
directed  electric  current  density  is  more  pronounced  for  the  disk  of  radius 
0.002A  than  for  the  disk  of  radius  0.02A. 

In  general,  all  the  E-field  solutions  in  [l]-[6]  begin  to  lose 
accuracy  somewhere  in  the  Rayleigh  region  and  become  worse  as  the  frequency 
decreases.  The  Rayleigh  region  [8]  is  the  range  of  frequencies  for  which 
the  maximum  dimension  of  S  is  much  smaller  than  the  wavelength.  These  E- 
field  solutions  fail  in  the  low  frequency  portion  of  the  Rayleigh  region 
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as  shown  by  the  following  reasoning.  When  the  frequency  is  sufficiently 
low,  the  magnetic  vector  potential  contributions  to  the  elements  of  the 
moment  matrix  are  insignificant  compared  with  the  electric  scalar  poten¬ 
tial  contributions.  As  a  result,  the  magnetic  vector  potential  contri¬ 
butions  are  lost.  The  remaining  scalar  potential  contributions  depend 

only  on  V  •  J.  Here,  V  •  J  is  the  surface  divergence  of  J.  Knowl- 
s  s 

edge  of  Vg  •  J  is  not  sufficient  to  determine  J.  Therefore,  the  E-field 
solutions  in  [l]-[6]  fail  in  the  low  frequency  portion  of  the  Rayleigh 
region. 

It  may  be  possible  to  obtain  reasonably  accurate  values  of  J  in  the 
low  frequency  portion  of  the  Rayleigh  region  in  the  following  manner. 

Usually,  there  is  at  least  one  frequency  that  is  high  enough  such  that 
one  of  the  E-field  solutions  for  .J  in  [ 1] — [ 6 ]  is  accurate  but  low  enough 
such  that  J  can  be  approximated  by  the  first  term  of  its  Rayleigh  series. 

The  Rayleigh  series  for  .J  is  its  low  frequency  expansion  in  non-negative 
integer  powers  of  the  frequency  [8].  If  there  is  such  a  frequency,  the 
first  term  of  the  Rayleigh  series  for  J  can  be  extracted  from  one  of  the 
E-field  solutions  in  [1]— [6] .  Knowledge  of  the  first  term  of  its  Rayleigh 
series  amounts  to  knowledge  of  J  in  the  low  frequency  portion  of  the  Rayleigh 
region. 

In  many  cases,  it  is  of  interest  to  obtain  not  only  J  but  also  the 
scattered  field.  The  scattered  field  is  the  electromagnetic  field  pro¬ 
duced  by  .J.  In  the  resonance  region,  knowledge  of  J  alone  is  sufficient  to 

calculate  the  scattered  field  because  the  density  q  of  electric  charge 

e 

associated  with  .J  can  be  calculated  from  the  equation  of  continuity 

V  •  J 
s 

q  =  : — 

e  -ju 


(1) 


In  (1),  0)  is  the  angular  frequency.  However,  in  the  low  frequency  portion 

of  the  Rayleigh  region,  knowledge  of  both  J_  and  qg  is  required  in  order  to 

calculate  the  scattered  field  because  qg  can  not  be  accurately  calculated 

from  .J  by  means  of  (1).  Accurate  calculation  of  qg  from  J_  by  means  of  (1) 

is  usually  not  possible  in  the  low  frequency  portion  of  the  Rayleigh  region 

as  shown  by  the  following  reasonirtg.  Usually,  J  approaches  a  solenoidal 

vector  function  as  the  frequency  approaches  zero.  As  a  result,  •  J_  is 

so  small  that  it  can  not  be  accurately  calculated  from  J.  Therefore, 

accurate  calculation  of  q  from  J  by  means  of  (1)  is  usually  not  possible 

e  — 

in  the  low  frequency  portion  of  the  Rayleigh  region.  However,  q^  can  be 

obtained  in  the  same  way  that  J  was  obtained  in  the  preceding  paragraph. 

In  the  manner  described  in  the  previous  two  paragraphs,  any  one  of 

the  E-field  solutions  [l]-f6]  for  J  can  be  extended  into  the  low  frequency 

portion  of  the  Rayleigh  region.  The  problem  of  the  perfectly  conducting 

surface  S  illuminated  by  an  incident  electromagnetic  field  in  the  low 

frequency  portion  of  the  Rayleigh  region  can  be  solved  by  other  methods. 

For  instance,  the  first  term  in  the  Rayleigh  series  for  J  can  be  obtained 

by  solving  a  magnetostatic  problem,  and  the  first  term  in  the  Rayleigh 

series  for  q  can  be  obtained  by  solving  an  electrostatic  problem.  In 
e 

most  of  the  literature  on  low  frequency  electromagnetic  scattering, 

Rayleigh  series  are  constructed  for  the  electric  and  magnetic  fields 
rather  than  for  the  electric  current  and  electric  charge.  For  example, 
see  [9]-[131. 

In  this  paper,  a  new  E-field  solution  is  presented  for  the  electric 
current  J  and  the  electric  charge  qg  induced  on  the  perfectly  conducting 
surface  S  immersed  in  an  incident  electromagnetic  field.  In  general,  S 


consists  of  several  surfaces  that  are  disjoint  from  each  other.  Some  of 
these  surfaces  may  be  open  and  others  may  be  closed.  The  new  E-field 
solution  is  similar  to  the  E-field  solutions  in  [1]- [6]  because  it  is  a 
moment  solution  of  the  electric  field  integral  equation  on  S.  However, 
the  new  E-field  solution  uses  different  expansion  and  testing  functions. 

In  the  new  E-field  solution,  the  expansion  functions  consist  of  two  sets 
of  vector  functions.  The  first  set  of  vector  functions  is  a  suitable 
basis  for  expanding  the  magnetostatic  current.  The  magnetostatic  current 
is  the  first  term  in  the  Rayleigh  series  for  J_.  The  second  set  of  vector 
functions  is  such  that  the  electric  charges  associated  with  them  form  a 
suitable  basis  for  expanding  the  electrostatic  charge.  The  electrostatic 
charge  is  the  first  term  in  the  Rayleigh  series  for  q^.  The  testing 
functions  in  the  new  E-field  solution  are  similar  to  the  expansion  func¬ 
tions.  These  expansion  and  testing  functions  render  the  moment  matrix 
well-behaved  as  the  frequency  approaches  zero. 

So  constructed,  the  new  E-field  solution  should  give  accurate  values 
of  J  and  qg  throughout  both  the  Rayleigh  region  and  the  resonance  region. 
Numerical  results  were  obtained  for  conducting  disks  of  radii  10  \  and 

0.021  and  for  conducting  spheres  of  radii  10  and  0.02X.  Each  disk  is 
excited  by  a  plane  wave  propagating  perpendicular  to  the  plane  of  the  disk. 
Each  sphere  is  excited  by  a  plane  wave.  These  numerical  results  agree  well 
with  the  known  solutions  for  the  electric  current  and  electric  charge  on  a 
small  disk  [14]  and  a  sphere  [15,  Eq.  (6-103)]. 
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II.  THE  MAGNETOSTATIC  CURRENT  AND  THE  ELECTROSTATIC  CHARGE 

Since  some  of  the  expansion  functions  for  the  new  E-field  solution 
will  form  a  suitable  basis  for  expanding  the  magnetostatic  current  and 
the  surface  divergences  of  the  rest  of  the  expansion  functions  will  form 
a  suitable  basis  for  expanding  the,  electrostatic  charge,  it  is  helpful  to 
establish  integral  equations  for  the  magnetostatic  current  and  the  electro¬ 
static  charge.  Such  integral  equations  can  be  obtained  by  substituting 

inc 

Rayleigh  series  for  ,J  and  I£  into  the  electric  field  integral  equation. 
The  electric  field  integral  equation  is  [4,  Eq.  (12)] 

-  -  E  (J)  =  -  |£nC  on  S  (2) 

H  —tan  —  n  —tan 

where  the  subscript  tan  denotes  the  component  tangent  to  S.  In  (2) , n  is 
the  impedance  of  free  space,  and  E(J)  is  the  electric  field  due  to  in 
free  space.  This  field  is  given  by  [16,  Sec.  2.1] 

»  1 

E(J)  -  -jn[k  JCr^GCr-r’Ms'  +  (V’  •  J(r '  ))G  (r-r '  )ds '  ]  (3) 

JJ  -  kjj  s - 


where 


G(r-r')  = 


3-jk|r-r» 
4tt  I  r-r '  I 


In  (3),  k  is  the  wave  number,  ds'  is  the  differential  element  of  area  at 

r'  on  S,  V'  •  is  the  surface  divergence  on  S  with  respect  to  the  coordi- 
s 

nates  of  _r' ,  and  r.  is  the  position  vector  of  the  point  at  which  E(J_)  is 

inc 

evaluated.  In  (2)  I£  is  the  incident  electric  field.  The  incident  elec¬ 
tric  field  is  the  electric  field  that  would  exist  if  S  were  absent. 

The  Rayleigh  series  for  J  is 


1  *  l  knJ 

n=0 


nT(n) 


where  J  is  independent  of  k.  It  is  assumed  that  the  surface  S  on 
which  J  resides  is  the  union  of  Q  surfaces  S^,  S^  which  are 

disjoint  from  each  other. 


s  ■  U  s, 

q=l  4 

In  (6) ,  \J  denotes  union.  The  surface  may  be  either  open  or  closed. 
If  is  closed,  then  it  has  no  edges.  If  is  open,  it  has  an  edge 
called  C  .  It  is  assumed  that  C  consists  of  R  closed  contours 

q  q  q 

{Cqr,.  r  =  Rq  ^  ' 


Rq 

C  =  O  C 

a  V  c 


Since  no  line  charge  can  accumulate  on  C  ^ ,  J_  must  satisfy 


J  •  u,  =  0  on  C 
—  — b  c 


q  -  1,2, ..  .Q 

r  «  1,2, .. . R 


where  u,  is  the  unit  vector  tangent  to  S  and  normal  to  C  .  It  follows 
-b  q  qr 

from  (5)  and  (8)  that 


u.  =  0  on  C 
— b  qr 


•  u^  =  0  on  C( 


The  ranges  of  values  of  q  and  r  in  (9)  and  (10)  are  the  same  as  in  (8). 

Similar  boundary  conditions  hold  for  J  ,  J  , ...  but  will  not  be  used 

inc 

The  Rayleigh  series  for  li  is 


ElnC  -  I  knE(n) 
n=0 


where  E'  is  independent  of  k.  The  expansion  of  G(r-r')  in  powers  of 
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DO 

<=(£-£')  ‘  £  I 


(~jk)n  | r-r ’ |n_1 


n=0 


n: 


(12) 


In  view  of  (5)  and  (12),  substitution  of  (3)  and  (11)  into  (2)  yields 


n-1  ,  ..n-m-1 


i  l"  I  S 


4ir  4  4  (n-m-1)! 

n=l  m=0 


Ls 


|r-r' |n  m'2  J(m)(r,)ds' 


tan 


°°  n+1  ,  ,  v  n-m+1 

+  J_  y  kn  y  lzl> _ y 

4tt  4  4  (n-m+l) !  Vs 

n=-l  m=0 


oo 

r-r ’ | n-m  y'  •  J(m)(r’)ds'  =  -  £  knE(n)(r) 


11  n=0 


(13) 


where  V  is  the  component  of  V  tangent  to  S.  The  operator  V  is  called 
s  s 

the  surface  gradient  on  S. 

Equation  (13)  implies  that  the  coefficient  of  kn  on  the  left-hand  side 
of  (13)  is  equal  to  the  coefficient  of  kn  on  the  right-hand  side  of  (13) 
for  n  =  -1,  0,  1,  2,...  .  Setting  to  zero  the  coefficient  of  1/k  on  the 
left-hand  side  of  (13),  we  obtain 


V'  •  J(0)(r') 

S  fY^H - 


ds '  =  0  on  S 


(14) 


In  view  of  (6) ,  (14)  implies  that 


J  J 
S 


V’  •  J(0)(r’) 

S  — 


r-r 


ds’  =  c(0q)  on  S  ,  q-l,2,...Q 

q 


(15) 


where  {c^q^}  are  unknown  constants.  We  set  S  =  S  ,  <J>  =  1,  and  W  = 


in  (A-l)  and  take  advantage  of  (9)  to  obtain 


V'  •  J(0)(r')ds'  =0,  q  =  1,  2,  ...Q  (16) 

JJ  s  - 

S 

q 

The  electric  charge  associated  with  j/^(r' )  by  means  of  (1)  is  called 


on  S 


8 


.(0) 


(*•')  = 


s  — 


J(0)(r') 


(17) 


-JO) 


Now,  (15)  states  that  each  of  the  surfaces  ( S  }  is  charged  to  a  constant 

q 

potential  whereas  (16)  states  that  there  is  no  net  charge  on  any  of  the 
surfaces  {S  }.  In  this  case,  it  is  well-known  from  potential  theory  that 
q^^(_r')  must  vanish  everywhere  on  S.  As  a  result,  the  magnetostatic,  cur¬ 
rent  satisfies 


7’  •  .J(n)  (r’  )  =  0  on  S 

s 


( 1 


Setting  the  constant  with  respect  to  k  on  one  side  of  (13)  equal 
to  the  constant  on  the  other  side  ol  (13),  we  obtain 


(1  ,  ,. 

1  f(  q  (r  } 

"  n-  ds'  =  F,  (r)  on  S 


4 tig  '  s  ;  j  r  r-r 

o  " 


-tan  — 


(19) 


where  gq  is  the  permittivity  of  free  space  and  q^  is  the  electric  charge 
associated  with  kJ^  by  means  of  (1). 


a)  •  i<n(r’) 

«(lV>  -  - 


(20) 


In  (20),  c  =  l/(pe  )  is  the  speed  of  light.  If 
o 

7  x  E(0)  =  0 


(21) 


then  there  is  a  scalar  function  <£^(r)  such  that 


E(0)(r)  =  -  V<j>(1)(r) 


(221 


and,  with  the  help  of  (6),  (19)  will  reduce  to 

7—^ —  !  f  --T - Vr-  ds '  =  -  4i(l)  (r)  -  C(lr°  on  S  ,  q=l,2, .  .  .Q  (23) 

4rr  I  1 r-r  -  q 

o  "  - 

where  are  unknown  constants. 
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The  following  reasoning  establishes  (21).  Substitution  of  the 

inC 

Rayleigh  series  (11)  for  E  and  the  Rayleigh  series 


HinC  -  l  kVn) 

n=0 

inc 

for  the  incident  magnetic  field  H  into  the  Maxwell  equation 


(24) 


v  X  Einc  =  -  jkn  Hlnc 


(25) 


gives 


l  k°V  x  E(n)  -  -  jr i  i  knH<"'1) 
n=0  n=l 


(26) 


Equation  (21)  is  a  consequence  of  (26). 

We  set  S*S  ,  <J>=1,  and  W  =  in  (A-l)  and  take  advantage  of 


(10)  to  obtain 


V’  •  J (1  (r  )ds  =  0  ,  q  =  1,2, ...Q 

JJ  a  -  - 

S 

q 


(27) 


In  view  of  (20),  (27)  implies  that 


j  q(1)(r,)ds'  =  0  ,  q  =  1,2,...  Q  (28) 

S 

q 

Now,  the  auxiliary  equation  (28)  atones  for  the  unknown  constants  {C^^}  in 
(23)  so  that  the  pair  of  equations  (23)  and  (28)  suffices  to  determine  q^\ 
It  is  evident  that  q^  is  the  electrostatic  charge.  Accordingly,  (23)  and 
(28)  a.e  called  the  electrostatic  equations. 

Intending  to  determine  the  magnetostatic  current  we  let  W 

be  a  differentiable  vector  function  tangent  to  {S  }  and  integrate  over  S 

q  q 

the  dot  product  of  W  with  the  terms  proportional  to  k  on  both  sides  of  (13) 


to  obtain 
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W(r)  •  (r)ds  + 


W(0  •  V<})(_r)ds  =  — 


W(r)  •  E(1)(r)ds,  q=l,2,...Q 

(29) 


In  (29), 


A(0)  (r)  = 


J<0)te'>  ds1 

4tt  I  r-r  *  I 


and,  thanks  to  (18), 


<P  = 


4tt 


V 


•  J(2)(r’) 


s  — 


r-r 


ds' 


(30) 


(31) 


Equation  (29)  will  reduce  to  an  equation  for  alone  if  W(r)  is  chosen 

such  that 

rr 

W(r)  •  7<|>(r)ds  =  0,  q-l,2,...Q  (32) 


According  to  (A-l) ,  (32)  will  be  satisfied  if 


V  •  W  =  0  on  S  ,  q=l,2,...Q 

s  —  q 


and 


W  •  u,  =  0  on  C  , 
—  — b  qr 


|  q- 1 , 2 , • • »Q 
r=l,2,...R 


(33) 


(34) 


If  (33)  is  true,  then,  according  to  (B-l),  there  is  a  scalar  function 
u(r_)  such  that 

W(r)  =  n  x  7  u(r)  on  S  ,  q=l,2,...Q 
- -  s  -  q 


(35) 


where  n  is  the  unit  vector  normal  to  S.  If  (34)  is  also  true,  then  u(r) 
must  satisfy 

j"q=l  >2, .  .  .Q 

u(r)  =  U  on  C  ,  <  (36) 

qr  qr  ) r=l , 2, . . .R 


where  (U  }  are  unknown  constants. 

qr 
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In  view  of  (32),  substitution  of  (35)  into  the  right-hand  side 
of  (29)  yields 

i  *  * 

J  W(r)  •  A(0)(r)ds  =  £  (n  x  V  u(r))  •  E(1) (r)ds,  q=l,2,...Q  (37) 

J }  ri  s 

S  S 

q  q 

Application  of  (C-l)  to  the  integral  on  the  right-hand  side  of  (37)  and 
subsequent  use  of  (36)  give 


of  in  (38)  by  the  right-hand  rule  stated  prior  to  (39). 

Equations  (35),  (C-l) ,  and  (36)  transform  the  left-hand  side 
of  (40)  to 

-  [f  u  n  •  V  x  A(0)ds  +  U  ff  A(0)  •  u.dJ. 

Jj  ~  -  r=i  qr  Jj  -  ~9- 

S  C 

q  fir 

Stokes'  theorem  [17,  Eq .  (42)  on  p.  489]  transforms  the  above  expression 


*  •  -7  *  -  Vq  f|  n  •  V  v  *<«)„» 


qr  J 1 


Therefore,  (40)  becomes 


[f  u (n  •  V  x  +  n  •  H(0))ds  =  P  U 


|  (n  •  7  x  A n  •  H(0))ds, 
'qr  q=l,2,...Q  0 


In  (41),  u  is  any  differentiable  scalar  function  that  reduces  to  the  arbi¬ 
trary  constant  U  on  the  contour  C  for  r  =  1,2,...R  .  Now,  (41)  can  be 
qr  qr  q 

valid  for  all  such  functions  u  onlv  if 


,(0)  „(0)  „  -on 

n  •  »  XA  -  n  •  H  on  S  ,  q  =  i  ,  2 ,  .  .  .  Q 

-  -  o 


f  , ..  A(0).  |f  ..(0). 

i  n  *  .  ^  A  ds  =  r  •  ts  ds, 
J  ’  ~  J 

S  S 

qr  qr 


rq=l  ,2,  .  .  .Q 


Since  A  ^  is  given  bv  (3h),  the  pair  of  equations  (42)  and  (43) 


helps  to  determine  a  .  Equations  19)  and  (18)  state  that  there  is  no 

,  ....  ,  .  r(0)  1 1„  /n't  / 1  Q  \  O.  U  - 


electric  change  associated  with  d 


Wc  believe  that  (9),  (18),  (42),  and 


(43)  uniquely  determine  the  magnetostatic  current  ,T 


Accordinglv,  these 


equations  are  called  the  magnetostatic  equations. 
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Equation  (9)  is  true  because  no  line  charge  can  exist.  Equation  (18) 
is  similar  to  [11,  Eq.  (1.80)].  Equation  (42)  is  a  statement  of  the  well- 
known  fact  that  the  normal  component  of  the  total  magnetic  field  is  zero  on 
a  conducting  surface  [11,  Eq.  (1.90)].  Equation  (43)  can  be  more  directly 
obtained  in  the  following  manner.  The  line  integrals  over  C  of  the  terms 
proportional  to  k  on  both  sides  of  (13)  are 


A(0)  •  u£  d£  +  V(j>  •  u£  d£  =  ^ 


E(1)  •  u£  d£  (44) 


where  Av  and  $  are  given  by  (30)  and  (31),  respectively.  Since  C^r  is  a  closed 
loop,  the  second  integral  on  the  left-hand  side  of  (44)  vanishes.  Appli¬ 
cation  of  Stokes'  theorem  [17,  Eq.  (42)  on  p.  489]  to  the  remaining  inte¬ 
grals  in  (44) ,  subsequent  use  of  (39) ,  and  multiplication  by  j  give  (43) . 

Equation  (43)  is  not  commonplace  because,  as  will  be  shown,  (43)  is 
necessary  only  if  S^  is  open  and  is  bounded  by  two  or  more  closed  contours 

(C  }.  If  S  is  not  open,  then  it  is  closed  and  therefore  not  bounded  by 
qr  q 

any  contour.  Consequently,  there  are  no  surfaces  S^  so  that  (43)  is 

absent.  If  S  is  open  and  is  bounded  by  one  closed  contour  C  ,,  then  S  , 

9  ql  ql 

is  identical  to  S  .  As  a  result,  (43)  is  redundant  because  it  can  be  obtained 

q 

by  integrating  (42)  over  S^.  If  is  open  and  is  bounded  by  closed  con¬ 
tours  (C  ,  r=l,2,...  R  }  where  R  >2,  then 

qr  q  q  - 


s  =  w  s 

q  rV:  qr 

As  a  result,  the  integral  of  (42)  over  is  the  sum  of  equations  (43) 
for  r«l,2,...  R  .  In  this  case,  (43)  must  be  enforced  for  R  -1  values 

q  q 

of  r.  At  first  glance,  there  appears  to  be  more  surface  area  on  the 


right-hand  side  of  (45)  than  on  the  left-hand  side.  However,  portions  of 
surface  area  with  oppositely  directed  normal  vectors  cancel  each  other  on 
the  right-hand  side  of  (45)  so  that  (45)  is  true. 
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III.  CONSTRUCTION  OF  THE  NEW  E-FIELD  SOLUTION 

The  new  E-field  solution  is  a  moment  solution  to  the  electric  field 
integral  equation  (2).  The  new  E-field  solution  is  constructed  by  expanding 
the  electric  current  J  as 


N  N 

T  cm  Tm  m  re  Te  .  Te 

"  j-1  J  A  J  1  "3 


(46) 


m  6 

where  J.  and  J.  are  vector  functions  that  are  on  S  and  are  tangent  to  S. 
-J  “I 

According  to  (6),  S  consists  of  the  surfaces  -( } .  For  convenience,  J™ 
is  chosen  to  be  non-trivial  only  on  for  q  equal  to  the  single  value 


m(j).  Similarly,  J_.  is  chosen  to  be  non-trivial  only  on  S 


e(j)  ’ 


Furthermore, 


{J™}  is  a  suitable  basis  for  expanding  the  magnetostatic  current  and 

e  (1) 

(V  •  J.}  is  a  suitable  basis  for  expanding  the  electrostatic  charge  q 

s  -J 

Neither  J™  nor  J?  is  allowed  to  depend  on  k.  The  magnitudes  of  J1?  and  J? 

-J  -J  -J  -J 

should  be  comparable  with  each  other.  Because  of  (8),  J1?  must  satisfy 


J ,  •  u,  =  0  on  C 
-I  b  qr 


ri-l,2,...N 

»  <q=m(.i) 


(47) 


l' 


r=l , 2, . . . R 


m(j ) 


and  Jj  must  satisfy 


J?  •  u,  =  0  on  C 
-j  -b  qr 


In  view  of  (18),  J?  must  also  satisfy 


1=1,2,. ..N 
q=e(j) 
r=l,2, . . ,R 


(48) 


e  (j  ) 


V'  •  j'J(r’)  =  0  on  S  , .  ,  j=l,2,...N 
s  — j  —  m(j)  m 


(49) 


m  G 

In  (46),  Ij  and  1^  are  unknown  coefficients  to  be  determined.  In 
general,  these  coefficients  will  depend  on  k.  The  scale  factors  {kp^ }  in 
(46)  are  for  later  convenience.  Here,  p  is  a  length  so  that  kp^  is 

dimensionless.  The  exact  value  of  p.  is  not  critical.  However,  the  order 

3 

of  magnitude  of  p^  should  be  that  of  a  dimension  of  S.  It  is  evident  from 
(46)  that  the  expansion  functions  for  J  are 


{Jm  ,  j=l,2,...  N  }  (50a) 

-j  m 

and 

(kPjJ®  ,  j=l,2,...  Ng}  (50b) 

In  analogy  with  the  above  expansion  functions,  testing  functions 

(W?  ,  i=l,2, . . .  N  }  (51a) 

—I  m 

and 

{kp.W*  ,  1=1,2, . . .  Ng}  (51b) 

are  introduced  on  S.  Both  W?  and  wf  are  tangent  to  S.  W™  is  chosen  to 

—i  — i  — i 

be  non-trivial  only  on  S  ,,,.  W.  is  chosen  to  be  non-trivial  only  on  S  .... 

m(i)  -i  e(i) 

Neither  W1^  nor  is  allowed  to  depend  on  k.  The  magnitudes  of  and 
should  be  comparable  with  each  other. 

The  testing  function  W™  is  chosen  to  satisfy 

W?  •  u.  =  0  on  C  ...  ,  i=l,2,...  N  (52) 

— i  — b  m(i) ,r  m 

0 

and  is  chosen  to  satisfy 

W?  •  u  =  6  on  C  ...  ,  i=l,2,...  N  (53) 

— l— b  e(i),r  e 

Moreover,  (V  •  W?)  should  be  a  suitable  set  of  testing  functions  for  the 
s  — i 

electrostatic  equation  (23).  Unfortunately,  (wT)  can  not  be  a  suitable 
set  of  testing  functions  for  the  magnetostatic  equation  (42)  because  (42) 
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is  a  scalar  equation.  However,  taking  a  cue  from  (49),  we  require  that 

Vs-^(r)-0  on  Sm(1),  1-1,2,...  »„  (54) 

where  Vg  •  is  the  surface  divergence  with  respect  to  the  coordinate  of  r^. 
Thanks  to  (54)  and  (B-l),  there  are  scalar  functions  {u_^}  such  that 

wj(r)  =  n  x  vg  u^r)  ,  i  -  1,  2,...  N  (55) 

Furthermore,  (W?)  are  chosen  such  that  if  an  expansion  for  in  terms 

of  {J™}  is  entered  into  (40)  by  means  of  (30),  then  enforcement  of  (40) 
for  W  successively  equal  to  each  member  of  {w1?}  determines  the  coefficients 
in  the  expansion  for  J  ^ . 

The  symmetric  product  of  two  vector  functions  is  defined  to  be  the 
integral  of  their  dot  product  over  S.  If  (46)  is  substituted  into  (3)  and 
if  (3)  is  subsequently  substituted  into  (2),  then  the  symmetric  products  of 
(2)  with  each  of  the  testing  functions  (51)  form  the  matrix  equation 

(56) 

In  (56),  I™  is  the  column  vector  of  the  coefficients  (i?)  in  (46),  and  le 
is  the  column  vector  of  the  coefficients  { I  ^  }  in  (46).  Also,  Zmm,  Zem,  Zme, 
and  Z  are  submatrices  whose  ijth  elements  are  given  by 

zTT  =  jk  ds  wT(r)  •  ds’  J?(r ’ )G(r-r ' )  (57) 

J  J  J  J  J 

Sm(i)  Sm(j) 

off  r  r 

zf?  =  jkZp,  ds  W®(r)  ds’  Jm(r ’ )G (r-r ’ ) 

ij  i  JJ  “  JJ  ~ 

Se(i)  Sm(j) 


(58) 


(' 


Z..  =  jk  Pj  ds  W\(r)  •  jj  ds  (r_'  )G(j:-r  ) 


ij 


Sm(i) 


Se(j) 


+  jp . 

1  j 
S 


ds  W?(0 


JJ 


ds' (V’ 
s 


JT (r'))C(r-r’) 


m(i) 


e  ( j ) 


zee.  -  ik^o.o. 

1.1  11  J 

s 


ds  W^(r) 


f  ds'  (r ' )G (r-r ' ) 


e(i) 


’e(j) 


ff 


+  jkp . p .  j  !  ds  W?(r )  ♦  7  i  ds' (V'  •  Je (r ' ) )G (r-r ' ) 
i  1  JJ  “  JJ  s  —j 


<  6 


e(i) 


e(j) 


->Tn 

Moreover,  V  and  V  are  column  vectors  whose  1th  elements  are  given  by 


m  1 


vi  =  n 


ff 


j  ^ 


W*  •  EinCds 


(6 


m(i) 


e  kpi 
V.  = 


i  n 


„e  „inc, 

•  E  ds 


(6 


Se(i) 


In  the  transition  from  (2)  to  (57)  -  (62),  it  was  permissible  to  omit  the 
subscript  tan  because  W1?  and  are  tangent  to  S. 

Thanks  to  (A-l) ,  (52),  and  (54),  (59)  reduces  to 


z“*  =  jkVj  |  ds  wj(0  *  | 


.1 


ds'  17 (r ' )G (r-r ' ) 


(6 


S 


m(i)  e(j) 

Tn  view  of  (53),  application  of  (A-l)  to  (60)  yields 


,ee 


'ij 


Jkpi°j  j 

s 


ds  W  (r ) 


ds'  1? (r ' )G (r-r ' ) 

Jj  "-I  - 


e(i) 


e(j) 


-  ikn.o. 


ds  (V 


We ( r ) ) 


ff 


ds'  (V'  •  J  .  (r ' )  )f.  (r-r 
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Substitution  of  (55)  into  (61)  and  subsequent  application  of 
(C-l)  give 

R  *  »  v 

^  ■  n  ||  ui(7  x  -lnC)  *  1  ds  +  n  \  ui(E±nc  *u£)dil  (65) 

*Vi)  Sn(i),r 

It  follows  from  (52)  and  (55)  that  u.  is  constant  on  C  .  Hence,  (65) 

1  rali;,r 

becomes 


.m 


n  J 


.  .  m(i)  r  . 

u  (V  x  E1  )  •  n  ds  +  ~  l  K.  I  Einc 

i  -  -  n  r=l  ir  J  “ 


•  u^d£  (66) 


m(i) 


"m(i),r 


where  K  is  the  value  of  u,  on  C  ,,,  .  Application  of  Stokes'  theorem 

ir  i  mtij  »r 

to  the  line  integral  in  (66)  gives 


R  . . . 

Vi  “  "  n  ||  Ui(V  X  -lnC)  *  JL  ds  +  £  "l  Kir  ||  S  *  ?  *  EinCds  (67) 
/i  r=l  _ 


m(i) 


m(i)  ,r 


Equation  (25)  reduces  (67)  to 


V®  =*  Jk 


u±(H 


inc 


n)ds  - 


Vi) 

jk  l  K 

r»l 


ir 


Vi) 


* 

H; 

. 


m(i),r 


inc 


•  n  ds 


(68) 


It  is  evident  from  (55)  and  (61)  that  addition  of  a  constant  to  u^ 
should  not  affect  the  value  of  v”.  Addition  of  a  constant  to  u^  does  not 
affect  the  value  (68)  of  V™  because,  according  to  (45), 


m 


Vi) 

(i)  "  ^  sm(i),r 


(69) 


If  some  linear  combination  of  the  {u^}  were  equal  to  a  constant,  then  the 
corresponding  linear  combination  of  the  {W?}  of  (55)  would  be  zero.  In  this 


case,  the  matrix  of  the  superscripted  Z's  on  the  left-hand  side  of  (56)  would 
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be  singular.  To  avoid  a  singular  matrix,  the  {u^}  should  be  chosen  so  that 
no  linear  combination  of  them  is  equal  to  a  constant. 

Construction  of  the  new  E-field  solution  is  now  complete.  The  new 
E-field  solution  for  .J  is  given  by  (46)  where  the  coefficients  { I™ }  and  {i^} 
are  the  elements  of  the  column  vectors  I™  and  Ie  that  satisfy  (56).  The 
ijth  elements  of  the  submatrices  Zmm,  ZGm,  Zme,  and  Zee  in  (56)  are  given  by 
(57),  (58),  (63),  and  (64),  respectively.  The  ith  elements  of  the  column 
vectors  v”1  and  Ve  in  (56)  are  given  by  (68)  and  (62),  respectively. 


IV.  LOW  FREQUENCY  BEHAVIOR  OF  THE  NEW  E-FIELD  SOLUTION 

If  the  wave  number  k  is  sufficiently  small,  the  elements  of  Z1™1”,  Zee, 

v”1,  and  V®  in  (56)  are  proportional  to  k  whereas  the  elements  of  Zem  and  Zme 

2 

are  proportional  to  k  .  As  a  result,  (56)  is  well  approximated  by  the  pair 
of  equations 


mrao  -+m 

M  1 

zeeo|e 


-*tno 

V 


4eo 

V 


(70) 

(71) 


whenever  k  is  sufficiently  small.  Here,  Z™110,  V™0,  ZEe°,  and  Veo  are  the 
low  frequency  limits  of  Z™”,  V™,  Zee,  and  Ve,  respectively.  From  (57), 
(68),  (64),  and  (62),  we  obtain 


z“°.Jk  ||  dsw^r).  ||  <!.'  (72) 

Sm(i)  Sm(j) 

£ 

Vi°  "  jk  {}  ui(5<0)  *Jl>ds  *  1  >Kir  }}  H(0)  *  B  ds  (73) 

Sm(i)  Sm(i),r 


where  E  and  H  have  been  extracted  from  the  right-hand  sides  of  (11) 
and  (24),  respectively. 

Later  in  this  Section,  it  is  shown  that  (70)  is  the  matrix  equation 
that  appears  in  a  moment  solution  for  the  magnetostatic  current  It 

is  also  shown  that  (71)  is  the  matrix  equation  that  appears  in  a  moment 
solution  for  the  electrostatic  charge  q^\  Presumably,  the  matrix  equations 
for  the  magnetostatic  current  and  the  electrostatic  charge  can  be  solved 
easily.  If  this  is  true,  then  the  matrix  equations  (70)  and  (71)  can  be 
solved  easily.  Hence,  the  matrix  equation  (56)  can  be  solved  easily  when  the 
frequency  is  low,  and  the  solution  will  tend  to  give  the  magnetostatic  cur¬ 


rent  and  the  electrostatic  charge. 


A  moment  solution  for  the  magnetostatic  current 


now 


con¬ 


structed.  The  matrix  equation  that  appears  in  this  solution  will  be  (70). 
In  view  of  (30),  (40)  is  an  equation  for  J .  Upon  substitution  of  j kV™ 
for  W  and  m(i)  for  q,  (40)  becomes 


where  u^  is  a  scalar  function  that  satisfies  (55)  and  is  the  value 


of  u.  on  C 

i  m(i),r 


If  the  expansion 


(7 


l  lm°Jm 


j=l 


J  -J 


is  inserted  into  (30)  and  if  (30)  is  substituted  into  (76),  then  equa¬ 
tions  (76)  will  form  the  matrix  equation 


^mmo  ^mo  _  +mo 


(7! 


Here,  1™°  is  the  column  vector  of  coefficients  {l?°}  in  (77).  In  (78),  the 


elements  of  zmmo  and  v™°  are  given  by  (72)  and  (73),  respectively.  As 


expected,  the  matrix  equation  (78)  for  I™0  is  the  same  as  the  matrix  equa¬ 


tion  (70)  for  I™. 


A  moment  solution  for  the  electrostatic  charge  q^"*  is  now  con¬ 


structed.  The  matrix  equation  that  appears  in  this  solution  will  be  (71). 
Substitution  of  the  expansion 

N 


q(1)(r')  =  j-  J  l!°  p.(V  •  Je(r')) 

1  —  r>  *-•  -l  i  e  — -i  — 

J-l 


(7< 


into  the  electrostatic  equation  (23)  and  integration  of  the  product  of 


(23)  with  -(kp./n)V  •  W.  over  S  ...  produces 

i  s  — i  e(i) 


N 


J  s 


ds(V  •  w; (r)) 

S  — X  — 


V'  *  JT(r') 


ds' 


4tt  r-r ' 


e(i) 


’e(j) 


kpt  rt 


.(1) 


<j>v  '(r)  vs  •  W^(r)  ds 


e(i) 


(1 , e(i ) ) 


ff 


V  •  W  (r)ds,  i  =1,2, .  .  .N  (8( 
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Thanks  to  (53)  and  the  divergence  theorem  [17,  Eq.  (42)  on  p.  503],  the 
second  integral  on  the  right-hand  side  of  (80)  vanishes.  Next,  (A-l)  is 
applied  to  the  first  integral  on  the  right-hand  side  of  (80),  and  then  (22) 
is  used  to  replace  the  resulting  V<f>^  by  -E^.  In  view  of  these  con¬ 
siderations,  equations  (80)  form  the  matrix  equation 


zeeo  -jeo  _  +eo 


(81) 


Here,  le°  is  the  column  vector  of  the  coefficients  {l?°}  in  (79).  In 

J 

(81),  the  elements  of  Z  and  V  are  given  by  (74)  and  (75),  respectively. 

->-00 

As  expected,  the  matrix  equation  (81)  for  I  is  the  same  as  the  matrix  equa¬ 
tion  (71)  for  ^e. 


V.  APPLICATION  TO  A  SURFACE  OF  REVOLUTION 

If  the  surface  S  is  a  surface  of  revolution,  suitable  expansion 
Tn  6 

functions  { }  and  {kp^J^}  for  the  new  E-field  solution  can  be  ic:  stn,..<ad 

d) 

by  taking  linear  combinations  of  the  expansion  functions  (J  }  and  {Jp^} 
defined  by  [5,  Eqs.  (2)  and  (3)] 


T,  (t) 


J  ,  «  u 
~nj  — t 


j=l,2,...P-2 
n=0,±l ,±2, . . . 


P,(t) 


j'*'  =  u 

~nj  — <|>  p. 


j=l,2, . . .P-1 


n=0,±l,±2, . . . 


Here,  t  is  the  arc  length  along  the  generating  curve  of  S,  and  <j>  is  the 
azimuthal  angle,  t  and  <p  are  orthogonal  coordinates  on  S.  ^  and  u^  are 
unit  vectors  in  the  t  and  <p  directions,  respectively.  Assuming  that 
tl,t2**’tP  are  Points  on  the  generating  curve,  Tj(t)  is  the  triangle  func- 


tion  which  begins  at  t j ,  peaks  at  tj+^,  dtl“  ends  at  t^^-  P  *-s  t“e  dis¬ 
tance  from  the  axis  about  which  the  generating  curve  is  rotated.  (t) 
is  the  unit  pulse  function  whose  domain  extends  from  t^  to  t^+^.  Pj  is 
the  value  of  p  at  the  center  of  the  domain  of  P^.(t). 

To  indicate  dependence  on  n,  the  expansion  functions  for  the  new 
E-field  solution  for  a  surface  of  revolution  are  called  {jm.}  and  (kp.J6 } 

j“nj 

instead  of  {J™}  and  {kp^J^}.  Similarly,  the  testing  functions  for  the  new 

E-field  solution  are  called  {W111  }  and  {kp.W6}  instead  of  (w?)  and  {kp  W?}. 

— ni  i— ni  — 1  i~ 1 

The  superscript  m  stands  for  magnetostatic,  and  the  superscript  e  stands 
for  electrostatic.  In  [5],  the  testing  functions  are  the  complex  conjugates 
of  the  expansion  functions.  Accordingly,  the  testing  functions  for  the  new 
E-field  solution  are  chosen  to  be  the  complex  conjugates  of  the  expansion 
functions  for  the  new  E-field  solution. 

wT  -  jm*  (84) 

— ni  —ni 


kp.  W.  -  kp,  J  . 
i  -ni  i  -ni 


Here,  *  denotes  complex  conjugate. 


Since  Jm.  and  kp.J6.  are  linear  combinations  of  (82)  and  (83),  Jm. 
— nj  j-nj  ~nj 


and  kp  j  J6  ^  are  proportional  to  e^n<^. 


It  can  be  shown  that  the  field  due  to 


any  electric  current  proportional  to  e^n^  is  also  proportional  to  e^n<^.  Hence, 

the  symmetric  products  of  the  fields  due  to  Jm.  and  kp.J6,  with  W™,  and  kp.W6. 

— nj  j-nj  —pi  i— pi 

are  zero  for  all  values  of  p  except  p=n.  As  a  result,  the  matrix  equation  (56) 
disintegrates  into  many  "smaller"  matrix  equations,  one  for  each  value  of n  in  (82) 


?me  -|m  ^m 

"n  n  n 

=  ,  n*0,±l,±2, . . . 
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„m 


According  to  (49),  the  surface  divergence  of  must  vanish. 


m 


This  means  that  can  not  have  any  electric  charge  associated  with 


.m 


it.  In  order  to  construct  J  ,  as  a  linear  combination  of  the  functions 

~nj 

in  (82)  and  (83),  we  have  to  know  the  electric  charges  associated  with 
these  functions.  The  surface  density  of  charge  associated  with 
is  called  q V  and  is  given  by  the  equation  of  continuity. 


V 


uj  -j(D  s  -nj 

Since  J1",  is  given  by  (82)  and  V  •  is  given  by  (B-3),  (87)  becomes 
— nj  s 


(87) 


P,(t) 


nJ 


(- 


-jojp  Aj 


_  _J±I 


(t) 


V 


-)  e 


jn<J> 


(88) 


where  A,  is  the  distance  from  t.  to  t.,,.  The  surface  density  qy,  of 
j  3  J+l  nj 

electric  charge  associated  with  is  given  by 


4> 

qnj 

=  ~  V 
-JU)  s 

•  A 

-nj 

(89) 

which  becomes 

q*. 

nj 

=  -=5-  P. 
wPjP  J 

<t).3n* 

(90) 

Noting  that  has 

no  electric  charge  associated  with  it. 

we 

choose 

and 

*4j  • 

j=l,2,...P-l 

(91) 

k% 

■4j  • 

j=l,2,...P-2 

(92) 

From  (84)  and  (85),  the  corresponding  testing  functions  are 

w1" 

W0i 

=  4i  ’ 

i=l ,2, . . .P-1 

(93) 

kpi4i 

-  jc 

^i  ’ 

i=l,2, . . .P-2 

(94) 

Thus,  for  n=0,  the  expansion  and  testing  functions  for  the  new  E-field 
solution  are  the  same  as  those  used  in  [5]. 


Equations  (92)  and  (82)  imply  that  is  proportional  to  1/k. 

6  0 
However,  J  .  is  not  allowed  to  depend  on  k  in  Section  III.  If  J  .  de- 
— nj  -nj 

pended  on  k,  it  might  be  difficult  to  obtain  an  accurate  numerical  solu¬ 
tion  to  (86)  because  (86)  would  not  be  properly  scaled.  On  the  other 
hand,  Z^m  and  Zq6  are  exactly  zero  so  that,  for  n=0,  (86)  separates  into 


nun  00 

two  matrix  equations,  one  involving  Z^  ,  the  other  involving  Z^  .  In 

this  case,  it  does  not  matter  how  Z^e  is  scaled  with  respect  to  Z^™. 

“*rn 

In  order  to  calculate  Vq  from  (68),  a  scalar  function  u^.,  must 


Oi 


be  found  such  that 


If 


Soi  '  a  *  Vs  u0i 


n  =  u,  x  u 
—  — <p  — t 


(95) 

(96) 


then  it  is  not  difficult  to  show  that  (95)  is  satisfied  by 


(97) 


Equation  (68)  was  derived  to  show  that  V™  is  proportional  to  k  for  small  k. 
inc 

If  E  is  such  that  the  integration  on  the  right-hand  side  of  (61)  yields 
an  expression  that  is  explicitly  proportional  to  k  for  small  k,  then  (68) 
is  not  necessary.  As  is  evident  from  [5,  Eq.  (80)],  such  is  the  case  for 
n»0  and  for  an  obliquely  incident  plane  wave.  Thus,  for  plane  wave  inci¬ 
dence,  the  matrix  equation  [5,  Eq.  (6)]  is  adequate  for  n=0. 
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where 


If  n  4  0,  it  is  evident  from  (88)  and  (90)  that 


t  .  <t>  d> 

q  .  +  a  .q  .  -  a  . . ,q  . . , 

nj  nj  nj  n.j+l’n.j+l 


nj  nA, 


Therefore,  it  is  suitable  to  choose 


Tm  TC  _L  7$ 

J  .  =  J  .  +  a  .J  . 
-nj  ~nj  nj-nj 


The  are  defined  by 


kp.J*.  =  kp.A  , 
j-nj  3~n3 


3=1,2, . . .P-2 


a  ...  Jy  ...  ,  < 
n,3+l  -m.3+1  ) n=±l ,±2 , . . . 


j=l,2, . . .P-1 


n=±l,±2, . . . 


(100) 


(101) 


From  (84)  and  (85),  the  corresponding  testing  functions  are 


tIm  Tt*  _<{>*  ,  T(J)* 

w,=J  -a,J.  +  a  . . ,  J  . . ,  , 

'“ni  — ni  ni— ni  n,i+l— n,i+l 


koi!£i  ■  kpi4i  ■ 


i=l,2, . . .P-1 


n=±l,±2, . . . 


In  view  of  (96),  it  is  not  difficult  to  show  that 


i=l,2,...P-2 
n=±l,±2, . . . 


(102) 


(103) 


where 


W1"  =  n  x  V  u  , 

-ni  —  s  ni 


"ni  '  n  Ti(t)  e'jn* 


(104) 


(105) 


Since  the  expansion  and  testing  functions  for  the  new  E-field 

solution  are  linear  combinations  of  the  expansion  and  testing  functions 

used  in  [5],  the  elements  of  the  superscripted  Z  's  in  (86)  are  linear 

n 

combinations  of  the  elements  of  the  superscripted  Z^'s  in  [5,  Eq,  (6)]. 
Of  course,  all  scalar  potential  contributions  to  the  elements  of  the 


superscripted  Z^'s  in  [5,  Eq.  (6)]  must  be  suppressed  from  the  calculation 

of  the  elements  of  Zmm,  Z6m,  and  Zme.  Otherwise,  severe  roundoff  error 

n  n  n 

would  occur  for  small  k.  Because  the  elements  of  the  superscripted  Z^’s 

2 

in  (86)  become  proportional  to  k  or  k  when  k  is  small,  machine  underflows 

will  occur  if  k  is  too  small.  In  an  attempt  to  avoid  underflows,  the  ele- 

1  2 

ments  of  the  superscripted  Z  's  in  (86)  were  divided  by  (x  kA  )  .  However, 

n  /i 

this  normalization  only  decreases  the  value  of  k  at  which  underflows  begin 
to  occur.  Unfortunately,  as  k  approaches  zero,  the  increasing  disparity 
between  the  real  and  imaginary  parts  of  the  matrix  elements  eventually 
exceeds  the  dynamic  range  of  the  computer. 

VI.  NUMERICAL  RESULTS 

The  new  E-field  solution  was  used  to  calculate  the  electric  cur¬ 
rent  J  and  electric  charge  qg  induced  by  a  plane  wave  axially  incident  on 
a  conducting  circular  disk  of  radius  0.002A  and  a  conducting  sphere  of 
radius  0.002A.  The  magnitudes  of  these  currents  and  charges  are  presented 
here.  The  new  E-field  solution  was  also  used  to  obtain  the  current  and 
charge  on  a  conducting  disk  of  radius  10  ^A  and  a  conducting  sphere  of 
radius  10  A.  These  currents  and  charges  are  not  shown  here.  It  suf¬ 
fices  to  state  that  they  agreed  well  with  the  known  currents  and  charges 
on  a  small  disk  [14]  and  a  sphere  [15,  Eq.  (6-103)].  The  new  E-field 
solution  could  not  be  calculated  for  disks  and  spheres  of  radii  consider¬ 
ably  less  than  10  ^A  because  of  machine  underflows. 

The  conducting  disk  lies  in  the  xy  plane.  On  the  disk,  the 
incident  field  (EinC,  HinC)  is  given  by 
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pine 

e  =  u  n 
—  —x 


(106) 


„inc 

H  =  -  u 

-  -y 


(107) 


where  u  and  u  are  the  unit  vectors  in  the  x  and  y  directions,  respec- 
-x  -y 

tively.  In  this  case,  can  be  expressed  as 


J  *  u  J  cos  <p  +  u, J,  sin  d> 
—  — t  t  — <p  <p 


(108) 


where  J  and  J,  are  dimensionless  functions  of  t.  In  (108),  t  is  the 
t  <p 

distance  from  the  center  of  the  disk,  <f>  is  the  azimuthal  angle,  u^  is 
the  unit  vector  in  the  radial  direction,  and  u,  is  the  unit  vector 
perpendicular  to  u^.  The  electric  charge  q^  can  be  expressed  as 

q  =  cos  <t 
e  c 


(109) 


where  c  is  the  speed  of  light  and  q  is  a  dimensionless  function  of  t. 

Figure  1  shows  |jt|  on  the  disk  of  radius  0.002X.  The  symbols 
x  represent  the  new  E-field  solution  for  | | .  These  symbols  are 
plotted  at  the  center  of  the  disk,  at  the  peaks  of  the  triangle  func¬ 
tions  {Tj  (t) }  in  (82),  and  at  the  rim  of  the  disk.  The  solid  curve  is 
the  known  solution  for  |j  |.  This  curve  was  obtained  by  calculating  the 
known  solution  at  31  points  equally  spaced  in  t  and  by  drawing  straight 
lines  between  them.  The  first  of  these  points  is  at  the  center  of  the 
disk.  The  last  point  is  at  the  rim  of  the  disk. 

Figure  2  shows  |j^|  on  the  disk.  The  symbols  x  represent  the 
new  E-field  solution  for  |j^|.  These  symbols  are  plotted  at  the  centers 
of  the  pulse  functions  {P^(t)}  in  (83).  The  solid  curve  is  the  known 
solution  for  |j^|.  This  curve  was  obtained  by  calculating  the  known 
solution  at  all  of  the  31  points  mentioned  in  the  previous  paragraph 
except  the  point  at  the  rim  and  by  drawing  straight  lines  between  them. 


Lng  circular 
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isk  of  radius 
the  new  E-field 
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The  known  solution  for  |j^|  approaches  infinity  as  the  reciprocal  of 
the  square  root  of  the  distance  from  the  rim.  Figure  3  shows  |q|  on  the 
disk  in  the  same  way  that  Fig.  2  shows  |j^|. 

In  Fig.  2,  the  x  nearest  the  rim  is  28%  higher  than  the  known 

value  of  |jjJ  at  the  corresponding  point.  In  Fig.  3,  the  *  nearest 

the  rim  is  36%  higher  than  the  known  value  of  |q|  at  the  corresponding 

point.  These  kinds  of  errors  are  to  be  expected  when  pulses  are  used 

to  expand  a  function  that  goes  to  infinity  as  the  reciprocal  of  the 

square  root  of  the  distance  from  an  edge  [18]. 

inc 

According  to  (107),  H  has  no  component  normal  to  the  disk. 
Therefore,  the  magnetostatic  current  vanishes  so  that  the  electric 

current  J[  reduces  to  kJ^  for  small  k.  In  Section  IV,  it  was  shown 
that  the  new  E-field  solution  for  J  approaches  as  k  approaches 

zero.  It  was  not  shown  that  if  vanishes,  then  the  new  E-field 

solution  will  reduce  to  kJ^  .  The  good  agreement  of  the  new  E-field 
solution  for  J  with  the  known  J  in  Figs.  1  and  2  is  fortunate.  It  is 
a  pleasure  to  state  that  the  new  E-field  solution  for  J_  agreed  just  as 
well  with  the  known  .J  on  the  disk  of  radius  10  *  A. 

The  conducting  sphere  is  placed  at  the  origin  and  is  illuminated 
by  the  incident  field 


_  inc  „  jkz 
E  =  u  n  eJ 
—  — x 


(110) 


„inc  jkz 

H  =  -  u  e 

-y 


(111) 


The  electric  current  and  electric  charge  qg  induced  on  the  sphere  are 
expressed  by  (108)  and  (109),  respectively.  For  the  sphere  of  radius 
0.002A,  Fig.  4  shows  |Jfc|,  Fig.  5  shows  |j^|,  and  Fig.  6  shows  jq|.  On 


INCIDENT  WAVE 
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an  the  conducting  sphere  of  radii 
Ls  x  represent  the  new  E-field  s< 
curve  is  the  known  solution. 


I q I  on  the  conducting  sphere  of  radius  0.002A.  Th 
symbols  x  represent  the  new  E-field  solution.  The 
solid  curve  is  the  known  solution. 


36 


this  sphere,  t  is  zero  at  (x=y=0,  z=-0.002A)  and  is  0.002ttA  at 

z=0.002A).  As  in  Figs.  1,  2,  and  3,  the  symbols  x  in  each 
of  Figs.  4,  5,  and  6  represent  the  new  E-field  solution  and  the  solid 
curve  is  the  known  solution. 


APPENDIX  A 


In  Appendix  A,  it  is  shown  that 


W 


V4>  ds  =  - 


4><V 


W)ds  +  j  'y(W  •  u^)d£ 
C 


(A- 


where  S  is  a  finite  surface,  ds  is  the  differential  element  of  area, 

V  •  is  the  surface  divergence  on  S,  4>  is  a  differentiable  scalar  func- 
s 

tion  defined  in  3-dimensional  space,  and  W  is  a  differentiable  vector 
function  defined  on  S.  Furthermore,  W  is  tangent  to  S.  The  surface  S 
may  be  either  open  or  closed.  If  S  is  closed,  then  the  second  term  on 
the  right-hand  side  of  (A-l)  is  to  be  omitted.  If  S  is  open,  then  S 
has  an  edge.  This  edge  consists  of  one  or  more  closed  contours  and  is 
called  C.  In  the  second  term  on  the  right-hand  side  of  (A-l),  d£  is 
the  differential  element  of  length  along  C  and  is  a  unit  vector  tan¬ 
gent  to  S  and  normal  to  C.  For  definiteness  u^  is  taken  to  point  away 
from  S. 


The  following  reasoning  is  used  to  show  that  (A-l)  is  true. 
Because  W  is  tangent  to  S,  the  divergence  theorem  [17,  Eq.  (42)  on 
p.  503]  for  <j>W  is 

\ 

0  ,  S  closed 


V  •  (<f>W)ds  = 
s  — 


(A- 


4>(W  •  u^)dg.  ,  S  open 


From  [17,  Eq.  (18)  on  p.  501], 


V  •  (AW)  =  Ar~  [A  (hnAW,)  +  A  (h,4>W0)l 


(A- 


where  (v^,  v^)  are  orthogonal  curvilinear  coordinates  on  S,  (h^,  V^)  are 
the  corresponding  metrical  coefficients,  and  (W^,  V^)  are,  respectively, 
the  components  of  W  in  the  directions  of  increasing  v^  and  v^.  Dif¬ 
ferentiating  the  products  on  the  right-hand  side  of  (A-3) ,  we  obtain 


Vg  •  (<j>W)  =  (J)Vg  •  W  +  W  •  Vg(J) 


(A-4) 


where 


VS  •  H  '  (3V-  <W  +  3^  (hlW2)! 


(A-5) 


HA  1  ,  1 

V  0  =  —  -7^-  u  +  —  •r-21'  U0 
s  h^  3v^  —1  ov^  —2 


(A-6) 


Here,  u^  and  u^  ara,  respectively,  the  unit  vectors  in  the  v^  and  v^ 
coordinate  directions.  According  to  [17,  Eq.  (18)  on  p.  501],  the  right- 
hand  side  of  (A-5)  is  indeed  the  surface  divergence  of  W.  From  [17, 

Eq.  (17)  on  p.  501],  the  right-hand  side  of  (A-6)  is  the  surface  gradient 
of  (j). 

Substitution  of  (A-4)  into  (A-2)  gives 


W  •  V  <J>  ds  =  - 
—  s 


<J>(V  •  W)ds  + 

s  - 


4>(W  •  u,  )dfc 
—  — b 


(A-7) 


S  SC 

Because  V  is  the  component  of  V  tangent  to  S  and  because  W  is 
s  ~ 

tangent  to  S, 


W  •  Vg(J)  =  W  •  V4>  on  S 


(A-8) 


Substitution  of  (A-8)  into  (A-7)  gives  the  desired  result  (A-l) . 
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APPENDIX  B 


In  Appendix  B,  it  is  shown  that  any  differentiable  vector  func¬ 
tion  W(r)  which  is  tangent  to  a  surface  S  and  which  has  no  surface 
divergence  can  be  written  as 

W(r_)  =  n  x  Vs<J>(r)  ,  £  on  S  ^B-l) 

where  n  is  the  unit  vector  normal  to  S,  Vg  is  the  surface  gradient 
on  S ,  and 


4>(r)  =  <j>(r  )  - 

— o 


(n'  x  W(r ’ ) )  •  dr' 


r  on  S 


(B-2) 


c 

In  (B-2),  r^  is  the  position  vector  of  an  arbitrary  point  on  S,  c  is 
any  contour  on  S  from  r^  to  £,  and  n'  is  n  evaluated  at  jc  * . 

The  following  reasoning  is  used  to  show  that  (B-l)  is  true. 
The  surface  divergence  of  W  is  called  Vg  •  W  and  is  defined  by  [17, 
Eq.  (18)  on  p.  501] 


\  >3^  <h2Wl>  +  3^  <hlW2)1 


(B-3) 


where  (v^,  v^)  are  orthogonal  curvilinear  coordinates  on  S,  and  (h^,  h^) 
are  the  corresponding  metrical  coefficients.  Also,  is  the  component 
of  W  in  the  direction  of  increasing  v^,  and  is  the  component  of  W  in 
the  direction  of  increasing  v^.  Since  W  has  no  surface  divergence, 

3^  (-h2Ul>  '  3^  <hlV 


In  view  of  (B-4),  the  differential  form 

hlW2dvl  ‘  h2Wldv2  (B_5) 

is  exact.  Therefore,  there  is  a  scalar  function  4>(v^,  V£)  such  that 


hlW2 


“Vi 


4>(vl,v2) 

4>(v1,v2) 


From  (B-6)  and  (B-7),  we  obtain 


(B-6) 

(B-7) 


W  = 


1  34  .  1  34 

~  h2  3v2  -1  h1  3vl  -2 


(B- 


where  (u^,  u2>  are,  respectively,  the  unit  vectors  in  the  directions 
of  increasing  v^  and  If  the  unit  vectors  (u^,  u2»  n)  form  a  right- 

handed  orthogonal  system,  then  (B-8)  can  be  rewritten  as 


where 


W  =  n  x  V  (p 
—  —  s 


V  <p  = 

3 


1  3d> 

-  - 1 —  11 

hL  dv±  -1 


+ 


1_  _3$_ 

h2  8v2  ~2 


(B- 


(B- 


According  to  [17,  Eq.  (17)  on  page  501],  the  right-hand  side  of  (B-10) 
is  indeed  the  surface  gradient  of  <f>.  Hence,  (B-9)  coincides  with  (B-l) 
It  is  evident  from  (B-8)  and  (B-10)  that 


V  (b  *  -  n  x  W 
s  —  — 


(B- 


The  desired  result  (B-2)  follows  from  (B-ll). 
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APPENDIX  C 


It  is  shown  in  Appendix  C  that 


(n  x  V  <J>) 


E  ds  =  - 


4> (V  x  E)  •  n  ds  + 


4>(E  *  <J  p  )  d  % 


(C-l) 


where  S  is  a  finite  surface,  ds  is  the  differential  element  of  area, 
n  is  the  unit  vector  normal  to  S,  is  the  surface  gradient  on  S,  <J> 
is  a  differentiable  scalar  function  defined  on  S,  and  E  is  a  differen¬ 
tiable  vector  function  defined  in  3-dimensional  space.  The  surface  S 
may  be  either  open  or  closed.  If  S  is  closed,  then  the  second  term  on 
the  right-hand  side  of  (C-l)  is  to  be  omitted.  If  S  is  open,  then  S  has 
an  edge.  This  edge  consists  of  one  or  more  closed  contours  and  is 
called  C.  In  the  second  term  on  the  right-hand  side  of  (C-l),  u^  is  the 
unit  vector  tangent  to  C  and  dH  is  the  differential  element  of  length 
along  C.  The  direction  of  n  is  the  direction  that  a  right-handed  screw 
would  advance  when  turned  in  the  direction  of  Uy  . 

The  following  reasoning  is  used  to  show  that  (C-l)  is  true. 
Stokes'  theorem  [17,  Eq .  (42)  on  p.  489]  for  <}>E  is 


0 


J  j  n  •  V  x  (<J>  E)  ds 
S 


<t>(E 


Uy)d£ 


S  c’osed 


S  open 


(C-2) 


If  (v^,  v^)  are  orthogonal  coordinates  on  S  and  if  (u^ ,  u^,  n)  form  a 
right-handed  system  where  (u^,  u2)  are,  respectively,  the  unit  vectors 
in  the  directions  of  increasing  v^  and  v2<  then  [17,  Eq.  (166)  on  p.  497] 


n  •  V  x  (<J>  E)  = 


hlh2 


k3v. 


(h24>E2)  - 


3v- 


(h1<()E1)) 


on  S 


(C-3) 


In  (C-3),  (h^,  h2)  are,  respectively,  the  metrical  coefficients  associ¬ 
ated  with  and  v2>  and  (E^,  E2)  are,  respectively,  the  components  of 
in  the  u^  and  u2  directions.  Differentiating  the  products  on  the  right 
hand  side  of  (C-3),  we  obtain 


n  •  V  x  (({>  e)  =  <|)(V  x  e^)  •  n  +  (n  x  Vg<J>)  •  on  S  (C-4) 


where 


1  3<*>  .  1  3<J> 

V  <p  =  : —  -  u.  +  i —  - —  u 


sy  h  3v  -1  h2  9v2  -2 


(C-5) 


According  to  [17,  Eq.  (17)  on  p.  501],  the  right-hand  side  of  (C-5)  is 


indeed  the  surface  gradient  of  <f> . 

Substitution  of  (C-4)  into  (C-2)  gives  the  desired  result  (C-l). 


43 


1 

! 

:  • 

» 

r  ' 

i 

: 

t 

5 

i 


REFERENCES 

[1]  A.  W.  Glisson  and  D.  R.  Wilton,  "Simple  and  Efficient  Numerical 
Methods  for  Problems  of  Electromagnetic  Radiation  and  Scattering 
from  Surfaces,"  IEEE  Trans.  Antennas  Propagat.,  vol.  AP-28, 

pp.  593-603,  Sept.  1980. 

[2]  M.  G.  Andreasen,  "Scattering  from  Bodies  of  Revolution,"  IEEE 
Trans.  Antennas  Propagat.,  vol.  AP-13,  pp.  303-310,  March  1965. 

[3]  J.  R.  Mautz  and  R.  F.  Harrington,  "Radiation  and  Scattering  from 
Bodies  of  Revolution,"  Appl.  Sci.  Res.,  vol.  20,  pp.  405-435, 

June  1969. 

[4]  J.  R.  Mautz  and  R.  F.  Harrington,  "H-Field,  E-Field,  and  Combined 
Field  Solutions  for  Conducting  Bodies  of  Revolution,"  AEU,  vol.  32, 
pp.  157-164,  April  1978. 

[5]  J.  R.  Mautz  and  R.  F.  Harrington,  "An  Improved  E-Field  Solution 
for  a  Conducting  Body  of  Revolution,"  AEU,  vol.  36,  pp.  198-206, 

May  1982. 

[6]  S.  M.  Rao,  D.  R.  Wilton,  and  A.  W.  Glisson,  "Electromagnetic 
Scattering  by  Surfaces  of  Arbitrary  Shape,"  IEEE  Trans .  Antennas 
Propagat . ,  vol.  AP-30,  pp.  409-418,  May  1982. 

[7]  C.  Cha  and  R.  F.  Harrington,  "Electromagnetic  Transmission 
through  an  Annular  Aperture  in  an  Infinite  Conducting  Screen," 

AEU,  vol.  35,  pp.  167-172,  April  1981. 

[8]  R.  E.  Kleinman,  "The  Rayleigh  Region,"  Proc .  IEEE  ,  vol.  53, 
pp.  848-856,  August  1965. 


[9]  A.  F.  Stevenson,  Solution  of  Electromagnetic  Scattering  Problems 
as  Power  Series  in  the  Ratio  (Dimension  of  Scatterer) /Wavelength, 


J.  Appl.  Phys.,  vol.  24,  pp.  1134-1142,  September  1953. 

[10]  R.  E.  Kleinman,  "Low  Frequency  Solution  of  Electromagnetic 
Scattering  Problems,"  in  Electromagnetic  Wave  Theory,  Delft 
Symposium,  pp.  891-905,  Pergamon  Press,  New  York,  1967. 

[11]  J.  Van  Bladel,  "Low-Frequency  Asymptotic  Techniques,"  in  Modern 
Topics  in  Electromagnetics  and  Antennas  by  H.  Bach  et  al., 

Peter  Peregrinus,  Ltd.,  Stevenage,  England,  1977. 

[12]  R.  E.  Kleinman,  "Low  Frequency  Electromagnetic  Scattering," 
in  Electromagnetic  Scattering  by  P.L.E.  Uslenghi  (ed.).  Academic 
Press,  New  York, 1978. 

[13]  T.B.A.  Senior,  "Low-frequency  Scattering  by  a  Perfectly 
Conducting  Body,"  Radio  Sci.,  vol.  17,  pp.  741-746,  July- 
August  1982. 

[14]  C.  J.  Bouwkamp,  "On  the  Diffraction  of  Electromagnetic  Waves  by 
Small  Circular  Disks  and  Holes,"  Philips  Res.  Rep.,  vol.  5, 
pp.  401-422,  December  1950. 

[15]  R.  F.  Harrington,  Time-Harmonic  Electromagnetic  Fields,  McGraw-Hill, 
New  York,  1961. 

[16]  C.  H.  Papas,  Theory  of  Electromagnetic  Wave  Propagation, 

McGraw-Hill,  New  York,  1965. 

[17]  J.  Van  Bladel,  Electromagnetic  Fields,  McGraw-Hill,  New  York,  1964. 

[18]  D.  R.  Wilton  and  S.  Govind,  "Incorporation  of  Edge  Conditions  in 
Moment  Method  Solutions,"  IEEE  Trans.  Antennas  Propagat.,  vol.  AP-25, 
pp.  845-850,  November  1977. 


